Deep Spitzer spectroscopy of the 'Flying Saucer' edge-on disk: 

Large grains beyond 50 AU. 
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>■ ! We present deep Spitzer-IRS low-resolution (A/AA ~100) 5-35 yum spec- 

: troscopy of the edge-on disk "the Flying Saucer" (2MASS J16281370-2431391) 

in the Ophiuchus molecular cloud. The spectral energy distribution exhibits the 
characteristic two-peak shape predicted for a circumstellar disk viewed very close 
to edge-on. The short-wavelength peak is entirely due to photons scattered off 
the surface of the disk, while the long- wavelength peak is due to thermal emis- 
sion from the disk itself. The Spitzer spectrum represents the first spectroscopic 
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detection of scattered light out to 15 yum from a bona-fide, isolated edge-on disk 
around a T Tauri star. The depth and the wavelength of the mid-infrared "val- 
ley" of the SED give direct constraints on the size distribution of large grains in 
the disk. Using a 2D continuum radiative transfer model, we find that a signifi- 
cant amount of 5-10 /zm-sized grains is required in the surface layers of the disk 
at radii of 50-300 AU. The detection of relatively large grains in the upper layers 
implies that vertical mixing is effective, since grain growth models predict the 
grains would otherwise settle deep in the disk on short time scales. Additionally, 
we tentatively detect the 9.66 /xm S(3) line of H2 and the 11.2 /xm emission feature 
due to PAHs. 

Subject headings: circumstellar matter - planetary systems: protoplanetary disks 
- stars: pre-main sequence - stars: individual (2MASS J16281370-2431391) - 
infrared: ISM 



Introduction 



Highly inclined circumstellar disks are particularly well-suited for studies of the dust 
properties in optically thick proto-planetary systems. The very high degree of obscuration 
toward the bright inner parts of these disks acts as a natural coronograph and thereby in- 
creases the contrast to the fainter outer parts by several orders of magnitude, enabling direct 
studies of the optical and infrared properties of the disk surface layers at radii of >50 AU; the 
region of the r = 1 surface to scattering as seen by an observer. In particular, the scattering 
properties and therefore the grain size distributions and the vertical height of the disk can 
be constrained using spatially resolved imaging. Examples of such s tudies employing opt i- 
cal and near-infr a red scattered light irnages o f edge-on disks include Burrows et al.l (Il996[) ; 
Stapelfeldt eraP Jl998h : ICotera et al.l J200lh : Istapelfeldt et all fl20o"3h : IWolf et airj2003h . 



These studies have shown that the extended emission at optical and near-infrared wave- 
lengths has a morphology closely resembling that theoretically expected for edge-on disks, 
i.e. two flattened reflection nebulosities bisected by a dark lane in the disk plane. Attempts 
to use these relatively short wavelengths to search for grains larger than those found in the 
interstellar medium have produced somewhat inconclusiv e results, and it has become clear 
that observations at longer wavelengths are required (e.g. iBrandner et al.ll2000l : IWood et al. 
2OO2I ) . A number of mid- infrared images have been made of edge-on di sks. Spatially resolve d 
emission at 11.8 /xm from the edge-on disk HK Tau B was presented by lMcCabe et al.l (120031 ) . 
who used the image to infer a grain size distribution with relatively large grains (1.5-3.0 /im). 
Perrin et al.l (120061 ) presented mid-infrared images of an edge-on disk clearly dominated by 
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Polycyclic Aromatic Hydrocarbon (PAH) emission rather than scattered hght. 

A problem is that edge-on disks tend to be very faint in the mid-infrared, requiring 
very sensitive observational facilities. Highly sensitive spectroscopy offered by the Spitzer 
Space Telescope provides new possibilities. Spitzer will not spatially resolve the mid-infrared 
emission, except for extremely large disks, but can obtain spectroscopic information that 
provides essential constraints for systems that have already been imaged at shorter wave- 
lengths. One important advantage of spectroscopy is that line emission can be distinguished 
from scattered light, a distinction that becomes essential when trying to constrain scattering 
properties of the dust. 

In this letter, we report deep mid-infrared spectroscopic observations using the Infrared 
Spectrograph (IRS) on Spitzer of the "Flying Saucer" edge-on disk located on the outskirts of 
the Ophiuchus star forming cloud core, first identified in high resolution near-infrared images 



(IGrosso et al.ll2003r). The spectrun i was obtained as part of the "From cores to disks" (c2d) 



Legacy program (lEvans et al.ll2003l ). The "Flying Saucer" is one of the most isolated edge-on 
disks known with a foreground visual extinction less than a few magnitudes, and is viewed 
almost exactly edge-on {i > 85°). 



2. Data reduction and results 

The "Flying Saucer" (2MASS J16281370-2431391, a =16^28™13!7, S =-24°31'39"[J2000]) 
was observed using both the short-low (SL2 = 5.2-7.5 /im, SLl = 7.5-14.5 /im) and long-low 
(LL2 = 14.2-21.0 /im, LLl = 19.8-38.1 /im) modules of the IRS (AOR KEY 0009829632). 
Each SL order was observed for 10 x 61s and each LL order for 20 x 31.5s. The spectral 
extraction was performed from the BCD pipeline (version S13.2.0) products by co-adding 
each nod position and then subtracting the two resulting 2D spectra to effectively remove 
the background. Due to the presence of a bright source on the peakup array for the SLl 
observation, the noise at wavelengths between 12.7 and 14.5 /im is significantly higher than 
in other parts of the spectrum. The "Flying Saucer" is clearly detected in SLl and 2 as 
well as in LLl; there is no detection in LL2, consistent with a location of the flux minimum 
around 20 /im. One-dimensional spectra were extracted using 3-pixel (5.4") apertures to 
optimize the signal-to-noise. To aid the modeling efforts, IRAC and MIPS photometry has 
been extracted from the c2d photometric database. The spectra match the IRAC band 3 
(5.7/im) and 4 (7.9 /im) photometry within the accuracy of the absolute calibration. The 
difference in slope between the spectrum and IRAC 3 and 4 may be due to the source being 
extended; the SL slit width is 3.6". The source is detected in MIPSl (24 /im) and only an 
upper limit is available in MIPS2 (70 /im) (see Tabled]). 
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Fig. 1. — The short-low spectrum showing the tentative emission features at 9.7 and 11.2 /im. 
The vertical bars show the location of common PAH features (solid bars) as well as the S(2), 
S(3) and S(4) lines of molecular hydrogen (dashed bars). The filled circles show the IRAC 
photometry in bands 3 and 4. The thick curve shows a spectrum that has been smoothed 
with a two-point box car filter. 
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The resulting SL spectrum is shown in Fig. [T] and shows two weak, narrow emission 
features at ~ 9.7yum and ~ 11.2 /im with hne-to-continuum ratios of ~ 1. Care should be 
taken in interpreting these features because the behaviour of the spectrograph is not yet 
fully understood, and features at the 5-10% level may be spurious. Assuming gaussian noise, 
the two emission features are detected with S/N ratios of 4-5. Inspection of the co-added 
2D spectra confirms that the features are compact and centered on the continuum position, 
at least along the slit, which is aligned at a position angle of 5°, or roughly perpendicular to 
the disk plane. We tentatively identify the emission features as the H2 S(3) line at 9.66 /im 
and the out-of-plane bending mode of PAHs at 11.2yum. The strength of the S(3) line 
corresponds to an integrated flux of (9 ± 2) x 10^^^ ergs~^ cm~^, while the 11.2 fim line flux 
is (15 ± 3) X 10~^^ergs~^ cm~^, using a 1st order polynomial fit to the continuum. For an 
excitation temperature of ~ 1000 K the other H2 lines are expected to be below the detection 
limit, taking into account the higher noise at the 6.9 //m S(5) line. The PAH line fiux is 2-3 
orders of magnitude less and the line-to-continuum ratio is a bout 10 times larg er than those 
determined for other PAH-emitting disks in the c2d sample (IGeers et al.ll2006l ). 

The spectral energy distribution (SED), shown in Fig. O declines between 4 and 14 /xm. 
The LL spectrum falls below the detection limit of ~ 1 mJy below 24 fim, but rises steeply at 



longer wavelengths. This spectral shape closely r esembles that t 



leoretically expected for a 



circumstellar disk viewed close to edge-on (see e.g. IWhitney et a. 



2003). The observed mid- 



infrared fluxes of 1 mJy matches closely those predicted by the iGrosso et al.l (120031 ) model. 
At short wavelengths (A < 40 /im) the optical depth to the source (the star and hot dust in 
the inner disk) is much greater than unity. This means all the short- wavelength photons are 
those that are scattered off the surface of the disk at large radii - in this case, 50-300 AU. The 
SED of the scattered light will have a shape that can be calculated roughly as the spectrum 
of the inner star-|-disk system multiplied by the albedo of the dust in the outer disk. At long 
wavelengths, the SED is dominated by thermal emission from the cold, outer part of the 
disk. These two components create a characteristic double-peaked SED with a minimum in 
the range 5-20 /xm that is unrelated to, but may in some cases be confused with, the silicate 
band at 9.7 /xm 



3. Constraining the dust size distribution 



We have constructed a 2D continuum radiative t ransfer model of the "Flying Saucer" 
with a setup similar to that of iPontoppidan et al.l (120051 ) using the Monte Carlo code RADMC 
( iDuUemond fc Dominikll2004l ). This setup uses an axisymmetric density structure on a polar 
grid to calculate images and spectra. We use Mie theory to calculate the input opacities to 
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Fig. 2. — Model fit to the SEP of th e "Flying Saucer" . The optical and near-infrared points 
are taken from iGrosso et al.l (120031 ). Grain size distributions weighted toward increasing 
grain sizes move the mid-infrared minimum in the SED to progressively longer wavelengths. 
The dashed curve the same as model 3, but with a proper implementation of anisotropic 
scattering. The dotted curve shows the SED with the scattered photons removed. No 
photons directly from the central star are seen at any wavelength. The colors of the curves 
refer to the same grain size distributions as those in Fig. [31 



Table 1: Mid- infrared photometry of the "Flying Saucer" 



Wavelength (/im) Flux [mJy] Reference 



3.6 


2.43 ±0.24 


c2d 


4.5 


2.32 ±0.23 


c2d 


5.7 


1.70 ±0.17 


c2d 


7.9 


0.84 ±0.08 


c2d 


24 


2.0 ±0.5 


c2d 


70 


< 300 


c2d 
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the model and explore a range of different grain size distributions. These are modeled by a 
power law size distribution dn{a)/da oc with maximum and minimum gr ain radii, amin 
andomax and a = 3.5. The dust material is silicate wi th optical co n stant s from lDraine &: Lee 
( 1984 ) and amorphous carbon (800 K experiment of Jager et al. ( 19981 )). For simplicity, it 
is assumed that carbon and sihcate grains are separate populations with 15% of the mass in 
carbon. Fig. [3] shows the scattering properties for some of the grain models. 



We model the physical disk structure using the iGrosso et al.l (120031 ) model parameters 
as a starting point. These are an outer disk scale height of if(100 AU)/100 AU = 0.15, 
a disk mass of 2 x 10~^ Mq, a surface density power law oc and a flaring scale 

height that varies as: h oc R^^^. The dis k has a radius of 300 AU. The disk parameters 



have been obtained by iGrosso et al.l (120031 ) by detailed fitting to the near-infrared images, 
using a grain size distribution with a maximum grain size of ~ 50 fim. Since no photometry 
is available above 70 pm, the properties of the disk mid-plane are not well-constrained. 

Some changes to the IGrosso et al.l (120031 ) m odel were mad e. First, we adopted a slightly 

smaller distance of 125 pc instead of 140 pc (Ide Geus et al.lll989l ). Second, to improve the 
fit to the near-infrared and IRAC photometry, we found it necessary to include a relative 
accretion luminosity of L^cc/L^ = 0.6 to the system, modeled by adding a blackbody with 
a temperature of 1400 K to th e input s tellar spectrum, which is modeled by a 0.084 
3500 K stellar atmosphere from iKuruca (119791 ) . as well as setting the foreground extinction 
to Ay = 0.5 mag. The accretion luminosity. Ay and effective temperature are somewhat 
interdependent, so other combinations may be possible, in particular, the central star may 
be cooler and more luminous, requiring less accretion. Removing the accretion component 
entirely will lower the 1-5 /im flux by a factor of ~ 2. In any case, these parameters do not 
affect wavelengths longer than ~ 5 /im, and therefore not the derived grain size distribution. 

Given a set of parameters determining the physical structure of the disk the grain size 
distribution strongly affects the shape and the location of the minimum of the SED in the 
mid-infrared. In Fig. [2], models with different grain size distributions are compared to the 
observed SED and in particular to the IRS spectrum. As larger grains are included in the 
grain population, the mid-infrared flux due to scattered photons increases relative to the rest 
of the SED. While light is always scattered at the surface where the scattering optical depth 
(^sca) is unity, the albedo controls how much light actually survives (is not absorbed and re- 
emitted at longer wavelengths) to be scattered at Tsca = 1- According to Fig. [3l a significant 
population of ~ 10 /im grains is required in order to have an albedo = Qsca/ (Qabs+Qsca) ^ 0.1 
at 20 fim. Additionally, as the small grains are removed, the thermal emission peak shifts to 
longer wavelengths. 



We ran a grid of models varying the minimum and maximum grain sizes, while keeping 
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Fig. 3. — Albedos and the angle-averaged phase functions or "asymmetry parameters" 
(< cosO > = g) for different grain size distributions. In principle, the albedo is a measure 
of the ability of a grain population to scatter light into the line of sight without removing it 
through absorption and re-emission at longer wavelengths. 



0.01 0.10 1.00 10.00 

Minimum grain radius [//m] 



Fig. 4. — Goodness-of-fit surface for the calculated grid of models. The contours have been 
normahzed, such that the best-fitting model has a goodness-of-fit of unity. The contours 
shown are in steps of 40%, the lowest having a value of 1.4. The hashed pattern indicates 
the area where no reasonable models can been calculated because amm > Omax • 
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the density structure constant, to minimize the quantity / {Fi^^^odcii^) — Fiy^ohs{ty)Ydh' ■ Figure 
m shows the goodness-of-fit surface. Inspection shows that good matches to the SED are 
found for grain size distributions with a^i^ > 0.5 /xm and Omax ^ 10 yum. While the exact 
hmits to the grain size distribution are not strictly constrained from the SED alone, the 
presented model clearly excludes grain size distributions with no grains larger than ~ 5/im 
or distributions that are dominated (in mass) by grains that are < 0.5 /xm. The upper 
cutoff for the grain size distribution may be significantly larger. The required presence of 
large grains is very stable to changes in the density structure of the disk and slope of the 
grain size distribution. Whether there is a lower cutoff to the grain size is less secure. For 
instance, a shallower grain size distribution was found to, at least in part, lessen the need 
for a lower cutoff of grain sizes. Additionally, a larger disk mass m ay be required sinc e 



our best-fitting J-band opacity is ~ 3 times lower than that used by iGrosso et al.l (120031 ). 
Models with higher disk masses lessen the need for a lower grain size cutoff by shifting 
the thermal peak in the SED to longer wavelengths, but do not affect the scattering peak. 
We are therefore confident that large grains are present in the Flying Saucer, while the 
evidence for a lack of small grains is somewhat more ambiguous. This result strengthens 
similar concl usions obtained for other edge-on disks using less constraining data at shorter 



wavelengths (IWood et al.l l2002l : iMcCabe et al.l l2003l ) and provides strong evidence for the 



presence of grain growth in the outer disk of the "Flying Saucer" . 

The model grid uses isotropic scattering, thereby requiring much fewer photon packages 
than if anisotropic scattering is used. While not a proper treatment of scattering of large 
grains, computing time is drastically cut this way. To ensure that the conclusions are not 
affected, we also calculated the SED of the best-fitting model using a proper anisotropic 
scattering model. Anisotropic scattering has the effect of making the image of the disk more 
compact around the central star. The effect on the SED is small as seen in in Fig. [21 it 
changes the flux by ~ 50% at wavelengths dominated by scattered photons. This difference is 
very small compared to the 2-3 orders of magnitude dependence on the grain size distribution 
in the mid-infrared. 



4. Discussion 

The IRS spectrum of the "Flying Saucer" has provided important constraints on the 
nature of disk material in the surface layers {H/R ~ 0.1) of the disk at radii larger than 
50 AU in several different ways. First, dust grains in this part of the disk are much larger 
than grains in the interstellar medium, clearly implying significant grain growth. Second, 
tentative emission features from PAHs and H2 reveal clues about emission excited by UV 
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photon processes in the outer disk. Had the disk been viewed face-on, these faint features 
would have disappeared in the glare from the warm inner parts of the disk; at 10 fim the 
"Flying Saucer" would be 200 times brighter if viewed face-on. 

Eviden ce from silicate emission features has already indicated that grains have grown to 
several /xm (IBouwman et al.ll200ll : Ivan Boekel et al.ll2005l ) at small disk radii, within ~ 10 A U 



(Ivan Boekel et al.l 120041 ). in Herbig Ae disks, as well as T Tauri stars (IKessler-Silacci et al. 



20061 ). In contrast, mid-infrared observations of scattered light such as those presented in 
this letter can reveal the presence of large grains in the surface layers much further away 
from the central star. While millimeter cor itinuum studies have also found eviden ce for 
large grains in the outer regions of disks (e.g. iNatta et al.ll2004l : iRodmann et al.ll2006l ). such 



observations probe only the mid-plane. In the case of the "Flying Saucer" , we have shown 
that significant grain growth has taken place at large distan c es fro m the star (up to several 
100 AU), consistent with the conclusions of iMcCabe et al.l (120031 ) in the case of HK Tau 
B, a similar edge -on disk. When comparing this re sult to recent models of grain growth 
and settling (e.g. iDuUemond fc Dominikll2004j . |2005| ). an apparent problem appears. These 
models predict that grain growth to sizes > 10 /im followed by rapid settling to the disk 
mid-plane takes place on very short timescales (10^ — 10^ years). This will effectively remove 
larger grains from the disk surface on a time scale that increases with radius and decreases 
with grain size. The new observational evidence presented here shows that large grains are 
present in the disk surface not only close to the star, as found in other studies, but also at 
much larger distances. This implies that in both these regimes there is an efficient mechanism 
for preventing grains in the ~ 10 fim range from settling. 

Adding to the puzzle is the indication of the presence of very small grains, such as 
PAHs, in the same general region of the disk as the large grains. Even if the detection of 
PAHs in the " Flying Saucer" is o nly tentative, some edge-on disks clearly show very strong 
PAH features (IPerrin et aPbooeh . Very small grains will not settle, so the question they 
pose is why they have not been removed by the grain coagulation producing the observed 
population of large grains, or whether there is a process to replenish them. 

Another possibility for explaining the presence of PAH emission is that the 11.2 /im 
feature may be light from the inner parts of the disk scattered on the outer surface. This 
requires that the line-to-continuum ratio of the PAH emission from the inner disk is the same 
as that observed in the "Flying Saucer" . Observations of face-on disks have shown that such 



strong line emission is rare, but not impossible (iGeers et al.ll2006l ). 



The tentative detection of the S(3) line of H2 is interesting. The S(3) line dominates 



the H2 line spectrurn 



'or ga s at a temperature of a few hundred Kelvin. Detailed modeling 



by iNomura fc Millarl (120051 ) has shown that, in the presence of a significant UV excess, the 
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S(3) line flux from a disk very similar to that of the "Flying Saucer" is expected to be 
~ 9 X 10~^^ erg cm~^, with most of the flux originating in the uppermost layers of the 
disk at radii of 10-100 AU. If confirmed, the presence of the line can therefore be interpreted 
as an indicator of a strongly enhanced UV field in the upper layers. 

Finally, the possible presence of high contrast emission lines in the spectrum of the "Fly- 
ing Saucer" suggests that cdgc-on disks may be excellent targets for looking for fluorescent 
tracers in the outer disk that would not have been seen if the disk had been viewed face-on. 
This illustrates the need for sensitive space-born high-resolution mid-infrared spectroscopy 
in the future offered by SOFIA and the James Webb Space Telescope. 
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